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In ultrarelativistic heavy ion collisions, the multiplicity of multi-strange baryons per participat- 
ing nucleon increases with centrality in a different fashion for different systems and energies. At 
RHIC, for copper-|-copper (CuCu) collisions the increase is much steeper than for gold-gold (AuAu) 
collisions. We show that this system size dependence is due to a core-corona effect: the relative 
importance of the corona as compared to the core (thermalized matter) contribution varies and the 
contribution of a corona nucleon to the multiplicity differs from that of a core nucleon. <j) mesons 
follow - as all hadrons - the same trend, but the difference between core and corona multiplicity is 
relatively small, and therefore the CuCu and AuAu results are quite similar. This simple geometrical 
explanation makes also a strong case in favor of the validity of Glauber geometry in the peripheral 
regions of ultrarelativistic heavy ion collisions, which is crucial for understanding the early evolution 
of the system. 
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Even before the first relativistic heavy ion beam has 
been delivered, the enhancement of the production of 
strange particles has been considered as a possible sig- 
nal for the existence of a plasma composed of quarks 
and gluons (QGP) This enhancement may occur if 
due to compression the chemical potential of the up and 
down quarks becomes that large that the system creates 
preferably strange quarks and antiquarks which materi- 
alize finally into strange hadrons. Since then strangeness 
enhancement is one of the hot topics in the analysis ul- 
trarelativistic heavy ion collisions. 

In the meantime heavy ion experiments have re- 
vealed that the multiplicity per participating nucleon of 
(multi)strange baryons is up to 20 times larger than that 
observed in pp collisions per participating proton at the 
same energy and that this enhancement is strongly cen- 
trality dependent [3| . The observed multiplicities in cen- 
tral collisions of heavy systems at RHIC and SPS energies 
can be well described assuming that the hadrons are in 
statistical equilibrium at a temperature close to the crit- 
ical temperature predicted by lattice gauge calculations 
01 0, S IB] • The chemical potential which is obtained by 
a fit to the data is small and therefore the relation to the 
originally predicted enhancement [l[ is not evident. 

Despite of many efforts, the centrality dependence 
of this experimentally observed enhancement has not 
yet found a generally accepted explanation. In (grand) 
canonical statistical models the enhancement for sym- 
metric systems is not dependent on the centrality. It 
has been advocated that the increase with centrality may 
be due to finite size effects (canonical suppression), but 
this gives only a sizable effect for very small volumes 
0. To agree with data one needs to assume that the 
effective volume depends on the number of participants 
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which is in contradiction to the observation at 



lower beam energies and which has not found a physical 
explanation yet. It has also been shown that the cen- 
trality dependence can by parameterized as a sum of two 
terms, one proportional to the number of participants, 
Npart, the other proportional to the number of binary 
collision, A^bm, @, however no physical interpretation of 
this dependence has been given. 

Picking up an older idea Q recently it has been pro- 
posed that the dependence of the multiplicity on the cen- 
trality has a geometrical origin [10,]: In this model one 
distinguishes between a dense volume area, referred to 
as core, and a low density peripheral region, referred to 
as corona. Only the core participates in a collective ex- 
pansion and produces particles in a statistical manner, 
whereas the corona region consists mainly of hadrons pro- 
duced in nucleon- nucleon collisions. Employing this idea 
in the framework of one of the most advanced models for 
the simulation of heavy ion reaction, EPOS 1.9, it has 
been shown that the difference in the centrality depen- 
dence between AuAu at RHIC and PbPb at the SPS can 
be explained in a natural way with this core-corona hy- 
pothesis. Later the authors of ref. Q advanced a similar 
idea. 

The parameters for this model have been determined 
to describe the AuAu data. Therefore, in order to vali- 
date this approach, data for other systems are necessary. 
The recent analysis of the RHIC CuCu data at the same 
beam energy provides such a data set: in the core-corona 
approach the particle multiplicities for different systems 
are determined by the relative size of the core and the 
corona. This is given by geometry. Therefore, the cen- 
trality dependence of the particle multiplicities for CuCu 
can be predicted without any new parameter. 

The purpose of this article is twofold: a) to demon- 
strate that these CuCu data are reproduced by Epos 1.9 
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and b) to show that the physics of the centraHty depen- 
dence of the different (multi)strange hadrons in different 
systems can be understood in a simple model based on 
the geometrically determined relative size of corona and 
core. 

Without any calculation some basic features of the 
CuCu data find their natural explanation in this core- 
corona approach. The observed identical enhancement 
(as compared to pp collisions) for central CuCu and 
AuAu collisions is due to the fact that central collisions 
are to almost 100% core (and the cores in CuCu and 
AuAu behave identically at the same beam energy). Con- 
sequently, the slope of the enhancement as a function 
of Npart is larger in CuCu. With decreasing centrahty 
the corona contributes more and more, up to iVpart = 2 
which is pure corona (= pp). In pp collisions the sup- 
pression (as compared to high energy e+e" strings) of 
strange baryons increases with the number of strange 
quarks which it contains Therefore the production 
of multistrange baryons is much more important in the 
core (thermal production) as compared to the corona 
(string decay). Dividing the central yield in heavy ion 
reactions by the yield in pp collisions the enhancement 
factor increases with the number of strange quarks con- 
tained in the baryon. This experimental observation has 
been baptized strangeness enhancement although, phys- 
ically speaking, it is a due to strangeness suppression in 
pp. 

This picture also explains naturally why strangeness 
enhancement is more pronounced at SPS as compared to 
RHIC. There are two reasons: First, at lower energies 
the corona contribution relative to the core contribution 
is larger over the whole range of iVpart- Second, as has 
already been discussed quite a while ago [ll|, the aver- 
age string mass in pp collisions is rather small and even 
at ^/s — 200 AGeV more than 30% of the strings have 
an energy smaller than the threshold for f2 production. 
Therefore the multi(strange) baryon/pion ratio is sup- 
pressed as compared to a very energetic e+e~ string, and 
this suppression increases with decreasing proton energy. 

In order to understand the physics of the centrality 
dependence of the multiplicities we develop first a toy 
model, which reproduces nevertheless quantitatively the 
experimental results. Later will we employ the sophis- 
ticated three-dimensional implementation of the core- 
corona concept in the EPOS model, which confirms the 
conclusion of this simple model. 

In our toy model, we separate the wounded (or partic- 
ipating) nucleons, Npart, into two classes. Those which 
have only scattered once are considered as "corona par- 
ticipants", which hadronize like a string in pp collisions. 
The others form the fireball and the particle multiplic- 
ity is given by phase space. Defining f{Ncoi-c) to be the 
ratio of participating "core nucleons" to all participants, 
and assuming that the multiplicity per wounded core nu- 
cleon, Mcoro: E^nd the one per wounded corona nucleon, 
-^corona, docs neither depend on Npart nor on the system 
size, the centrality dependence of hadron multiplicities is 



given as 

(iVpart) (1) 
= ^part [/(iVcorc) ' M^„,, + (1 - /(iVcorc)) ' M^„„ J 

where i refers to the hadron species. It is important to 
notice that for a given total number of wounded nucle- 
ons, / depends on the size of the interacting system. For 
RHIC experiments at = 200 and 62 AGeV, /(A^core) 
has recently been published [l3]. However, there the cal- 
culations are done for the center of the centrality bins, 
whereas here we use a Glauber model calculation which 
averages f{Ncorc) over the centrality bins. This gives dif- 
ferent results in more peripheral reactions. Mcoro can ei- 
ther by taken from statistical model fits or can be directly 
determined from experiment. For our calculation, we fix 
Afcorc by applying eq. [T]to the most central Au-|-Au or 
Pb-I-Pb data point, by replacing the left hand side of eq. 
[1] by the experimental value. Mcorona is given as half of 
the multiplicity measured in pp collisions. Once these 
parameters are fixed, the centrality dependence of M* is 
determined by eq. [1] Especially the centrality depen- 
dence of the lighter Cu-|-Cu system follows then without 
any further input. 
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TABLE L Corona, Mcorona, and bulk, Mcoro, multiplicity per 
wounded nucleon as well as statistical model predictions 0] 

for A/corc 

The values of Mcoro and Mcorona as well as the thermal 
predictions Q of Mc^* are summarized in table 1. All 
multiplicities in this paper refer to dn/dy\y^o. To deter- 
mine Mcorona we use the data of [l^ for = 200 AGeV 
and EPOS 1.9 results for = 62 AGeV. Mcorona is de- 
termined from pBe for the NA57 and the NA49 data. The 
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experimental results can be found in [13, llJ, ll5|j UM ilZl 
and are partially still preliminary. We include in our 
analysis the centrality dependence of the (j) meson be- 
cause in the past it has been argued that it reflects the 
increased production of s-quarks in a dense medium . 

The results of our calculation (eqUJ) are displayed in 
fig. 1 for reactions at = 200 AGeV, = 62 AGeV, 
and -^/i = 17.2 AGeV, respectively. We plot the so- 
called nuclear enhancement, i.e. the ratio of the multi- 
plicity per participant to half the multiplicity in proton- 
proton scattering, as a function of Nyjound ^ for differ- 
ent hadrons and different systems, ft is evident that the 
simple toy model, i.e. eq. [TJ describes the experimentally 
observed centrality dependence for all three energies, for 
the strange baryons as well as for the (f>. Of course the 
assumptions of an abrupt transition between string frag- 
mentation and fireball formation is rather crude but the 
large experimental error-bars make it impossible to re- 
fine this model. At SPS energies data sets from 2 ex- 
periments, NA49 and NA57, are available. They do not 
agree. It is therefore useful to display calculations for 
both data sets. Our calculation reproduces the NA49 
data but is for peripheral reactions less steep than the 
NA57 data set. 

It is interesting to display the core corona effect in a 
different way. In fig[5]we display the multiplicity per par- 
ticipating nucleon of all observed hadrons in peripheral 
reactions normalized to the same multiplicity in the most 
central collisions as a function of M^oronaf-^^core- We see 
that all, non-strange as well as strange, particles fall on 
the same curve. More precisely, both theory and experi- 
ment show exactly the same linear increase. This means 
that the enhancement of the multiplicity in central heavy 
ion collisions as compared to peripheral reactions is for 
all particles strongly correlated with the multiplicity dif- 
ference in pp collisions as compared to what is expected 
in a purely thermal environment. 

Let us now turn to a somewhat more sophisticated im- 
plementation of the core-corona picture in the framework 
of the EPOS approach [lO|. EPOS is a parton model, 
so in case of a nucleus-nucleus collision, there are many 
binary interactions, each one represented by a parton 
ladder [l9l[20|. Such a ladder may be considered as a lon- 
gitudinal color field, conveniently treated as a relativistic 
string. In case of proton-proton scattering, the strings de- 
cay via the production of quark-antiquark pairs, creating 
in this way string fragments - which are usually identified 
with hadrons. When it comes to heavy ion collisions, the 
procedure is modified: one considers the situation at an 
early proper time tq, long before the hadrons are formed: 
one distinguishes between string segments in dense areas 
(more than some critical density po segments per unit 
volume), from those in low density areas. The high den- 
sity areas are referred to as core, the low density areas as 
corona. Compared to our toy model, here we perform a 
real three-dimensional core/corona separation: the core 
contribution will be more important at central rapidities 
as compared to projectile and target rapidities. 
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FIG. 1: Enhancement of nuclear strange particle production 
for different energies/systems. We compare toy model pre- 
dictions (lines (dashed lines for the NA49 results)) and data 
(points). 



Whereas corona particle production is as in proton- 
proton scattering, the core part is assumed to be collec- 
tively expanding: Bjorken-like in longitudinal direction 
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FIG. 2: Ratio of the peripheral and central multiplicity of 
all particles observed in 200 AGeV Au+Au reactions as a 
function of the ratio of the corona and the core multiplicity. 



with in addition some transverse expansion. We assume 
particles to freeze out at some given energy density epo j 
having acquired at that moment a collective radial flow. 
The latter one is characterized by a linear radial rapidity 
profile from inside to outside with maximal radial rapid- 
ity J/rad (model parameter). Hadronization then occurs 
according to covariant microcanonical phase space. So 
the core definition and its hadronization are parameter- 
ized in terms of few global parameters, for details see 

m- 

In fig. 3, we plot the nuclear enhancement obtained 
from EPOS 1.9. We observe, as already in the simple toy 
model, a steeper increase of multiplicities with centrality 
in case of CuCu compared to AuAu, due to a different 
-^part dependence of the relative core/corona weight in 
the two systems. 

In conclusion, we have shown that the centrality de- 
pendence of strange hadron production can be well ex- 
plained as a purely geometrical effect: the total volume 
is composed of a low density part ~ the corona, and a 
high density part - the core. Particle production in the 
corona is essentially string breaking as in pp collisions, 
whereas the core represents a fireball from which hadrons 
are produced according to phase space as in statistical 
models. This approach explains in a natural way the 
striking difference in the centrality dependence in CuCu 
compared to AuAu, being simply due to a different rel- 
ative core/corona weight at a given A'part for the two 
systems. 
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FIG. 3: Enhancement of nuclear strange baryon production 
for different energies/systems. We show the full EPOS 1.9 
results as well as the core contribution (dotted). Data are 
displayed as points. 



The (f) mesons fit perfectly into this picture, the differ- 
ence with respect to other hadrons like or S is quanti- 
tative. The ratio of core to corona multiplicity is consid- 
erably bigger for the latter ones, therefore the enhance- 
ment increases faster with iVpart than in case of 0. As a 
consequence, the enhancement for CuCu is close to the 
AuAu result. 

We have demonstrated that a very simple core/corona 
toy model already explains the main features seen in 



5 



the data. More sophisticated EPOS simulations confirm 
these finding ^21j . 

It is a consequence of this observation that the produc- 
tion of strange baryons - as that of all other hadrons - fol- 
lows very simple rules: after subtracting the corona part 
(whose importance depends on centrality / system size), 
all particles are created from some bulk matter whose 
properties are independent of the size of the system. It 
is actually remarkable that the core part is strictly de- 
termined by geometry and becomes therefore very small 
for the very peripheral reactions. This continuity implies 
that the expected new state of matter would exists also in 



very small systems. Another important consequence of 
this observation is that the hadron yield of hadrons con- 
taining u, d, and s quarks only, is sensitive to the freeze 
out of the matter. If one wants to have information on 
earlier stages other probes have to be used. 
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